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COMPUTATION OF TONE NOISE FROM SUPERSONIC JET IMPINGING ON FLAT 

PLATES 


Ching Y. Loh 
Taitech, Inc. 

Beavercreek, Ohio 45430 


Abstract 

A supersonic jet impinging normally on a flat 
plate has both practical importance and theo- 
retical interests. The physical phenomenon is 
not fully understood yet. Research concen- 
trates either on the hydrodynamics (e.g. lift 
loss for STOVL) or on the aeroacoustic load- 
ing. 

In this paper, a finite volume scheme-the 
space-time conservation element and solution 
element (CE/SE) method-is employed to nu- 
merically study the near-field noise of an un- 
derexpanded supersonic jet from a converg- 
ing nozzle impinging normally on a flat plate. 

The numerical approach is of the MILES type 
(monotonically integrated large eddy simula- 
tion). The computed results compare favorably 
with the experimental findings. 

1 Introduction 

High speed impinging jets are important to aircraft in- 
dustry. For vertical landing and short take-off (STOVL) 
airborned vehicles propelled by jet engines, there are un- 
desirable and adverse effects from the impinging jets. 
These include up to 60% lift loss and the acoustic load- 
ing and noises generated by a feedback loop between the 
jet and the ground. 

Research on jet impinging normally on a flat plate has 
unfolded in a broad way. Some researchers concentrate 
on the hydrodynamics and investigate the lift loss, other 
researchers focus on the sound produced and the aeroa- 
coustic loading of the impinging jets. There are a vast 
number of papers published on these topics. For exam- 
ple, after careful observations over the experimental data, 
Powell [1] pointed out that the small instability waves 
(vortices) around the jet shear layers and the consequent 
radial wall jet are responsible for the noise as they inter- 
act with the flat plate and produce sound waves. Ho and 
Nosseir [2] explained the feedback loop in the imping- 
ing jets; while Tam and Ahuja [3] put forward another 


theoretical model for the acoustic feedback loop. In the 
meantime, experimentalists are conducting various phys- 
ical experiments in search for the feedback mechanism 
in different impinging jet situations. Alvi and Iyer [9] 
studied the impinging jet with lift plate, Krothapalli et al 
[8] investigate both the jets with convergent nozzle and 
convergent-divergent (C-D) nozzles. Henderson etal\ 4 - 
7] performed experiments of sound producing impinging 
jets on small and large plates. 

Numerical simulations of unsteady impinging jets 
have also been carried out. Kim and Park [10] used the 
popular TVD upwind scheme. Sakakibara and Iwamoto 
[11] also used TVD scheme to study oscillations in im- 
pinging jets and the generation of acoustic waves. 

In the present paper, a MILES (monotonically inte- 
grated large eddy simulation) type scheme is used to in- 
vestigate the problem of a jet impinging on a flat plate. 
The MILES approach appears somewhat similar to LES, 
but there is no explicit filtering since the cell- averaging 
process is already equivalent to spatial filtering. Due to 
certain built-in numerical dissipation in a MILES finite 
volume scheme, the SGS (subgrid scale) model is not 
needed. The recent space-time conservation element and 
solution element method (CE/SE) [12-13] is a MILES 
type finite volume method with generally less numeri- 
cal dissipation and is adopted for the computation. As 
demonstrated in previous papers, the CE/SE scheme is 
well suited for aeroacoustics computation [15,16]. Be- 
cause of the CE/SE non-reflecting boundary conditions 
(NRBC), which are based on the physics of plane wave 
propagation [14], a smaller near field computational do- 
main can be used in the present numerical simulation and 
helps to save both memory and CPU time. 

The governing equations and the 2-D axisymmetric 
unstructured Navier-Stokes (N-S) CE/SE scheme used 
here is briefly reviewed in Section 2. Section 3 illustrates 
the noise problems of an impinging jet on normal flat 
plates, both large and small, with the initial and bound- 
ary conditions. The numerical results are presented and 
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compared to available experimental findings [4 - 7] in 
Section 4. Concluding remarks are drawn in Section 5. 

2 The Governing Equations and the 
Unstructured Axisymmetric CE/SE 
Navier-Stokes Solver 

As our main concern lies in the aeroacoustical behav- 
iors of the impinging jets which remain axisymmetric in 
the problems under investigation according to the exper- 
imental results [4-7], it is appropriate to adopt and solve 
the axisymmetric Navier-Stokes equation system. 


Qz = -Ui/Uiy, Q 4 = -Gi/y. 

By considering (x,y,t) as coordinates of a three- 
dimensional Euclidean space, and using Gauss’ di- 
vergence theorem, it follows that Eq. (1) is equivalent to 
the following integral conservation law: 




m = 1, 2, 3, 4, (2) 


where S(V) denotes the surface around a volume V in 

. 2/3 and H rn (F m , Grm • 


2.1 Conservation Form of the Unsteady 
Axisymmetric Navier-Stokes Equations 

Consider a dimensionless conservation form of the un- 
steady axisymmetric Navier-Stokes equations of a per- 
fect gas. Let p,u,v,p, and 7 be the density, streamwise 
velocity component, radial velocity component, static 
pressure, and constant specific heat ratio, respectively. 
The axisymmetric Navier-Stokes equations then can be 
written in the following vector form: 


Ut + F x + G y — Q, (1) 


where x, y > 0 , and t are the streamwise and radial co- 
ordinates and time, respectively. The conservative flow 
variable vector U and the flux vectors in the streamwise 
and radial directions, F and G , are given by: 


U = 



F = 


f Fl \ 

f 2 

f 3 

\fJ 


G = 



with 

Ui = p, U 2 = pu, U 3 = pv, 

Ui = p /( 7 - 1 ) + p{u 2 + v 2 )/2. 

The flux vectors are further split into inviscid and viscous 
fluxes: 

F = Fi — F v , G = Gi - G v , 

where the subscripts i and v denote ‘inviscid’ and ‘vis- 
cous’ respectively. Details of these terms can be found in 
e.g. [17, 16]. 

The right hand source term Q is the same as in the 
axisymmetric Euler equations [17, 16]: 


2.2 Treatment of the Source Term 

The treatment is identical to the one used in [17] and is 
briefly reiterated here. Since the source term Q itself is a 
function of the unknown U , a local iterative procedure is 
needed to determine U. The discretized integral equation 
( 2 ) reduces to the form 

U -Q(U)At = U H , (3) 

where U h is the local homogeneous solution (Q = 0 
locally). Note that U h only depends on the solution at 
the previous time step, i.e. , U h is obtained using explicit 
formulas. A Newton iterative procedure to determine U 
is then 


U i+ 1) = £/ (i) - (!jj) -1 [*(U (<) ) - U H ], 


where i is the iteration number and 


*(E7) = U - Q(U)At. 


Normally, U at the previous time step is a good initial 
guess U ^ and the procedure takes about 2-3 iterations 
to converge. The Jacobian matrix is given by 


dU 



1 

0 

At 

y 

0 


U 2 U*At 

U(y 

1 | U3 At 

1 ^ u iy 

U 2 At 

u iy 

0 


[/.? At 

ulv 

0 

l U3 At 

i+ U x y 

0 


A 1 

A2 

a 3 

-j , jU 3 At 
1 ^ 2 U iy 


where 


u 3 (i-l)((/| + u|) 

1 = - ~u^ bv> u 


Q = 


Q2 1 
Qz 
\qJ 


where 


Q 1 — —Us/y, Q 2 — —U2Us/Uiy, 


As 


A2 

At 

y 


-(7-1) 


U 2 U 3 At 

U ly 


7 — 1 ^ 1+3 Uj jUa 

2 Ul U\ 


The inverse of the Jacobian, i.e., j ) -1 can easily be 

derived analytically for this particular case, thus, leading 
to a savings in CPU time. 
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2.3 Review of the CE/SE Numerical Scheme 

The space-time conservation element and solution ele- 
ment(CE/SE) method is a recently developed finite vol- 
ume method, with second order accuracy in both space 
and time. Here, CE (conservation element) stands for 
a control volume or cell, while SE (solution element) 
stands for the cell interfaces. Despite its nominal sec- 
ond order accuray, the scheme may offer high resolution, 
lower dissipation and dispersion errors [12, 13]. As a re- 
sult of the following advantageous features, the CE/SE 
scheme is chosen as the numerical tool for computation: 

1. conservation in both space and time, the integral 
equations of conservation laws are literally solved; 



Figure 1: A sketch of the impinging jet problem. 


2. only a compact cell stencil is needed, (hence both 
the conservative variables U and their gradients 
U x , U y are unknowns); 

3. careful and accurate surface flux calculation; 

4. truly multi-dimensional, simple but effective non- 
reflecting boundary conditions (NRBC); 

5. effortless implementation of computation, no nu- 
merical fix or parameter choice is needed; 

6. the scheme is of MILES type (monotonically in- 
tegrated large eddy simulation), the finite volume 
cell-averaging plays a role of filter, while the van 
Albada limiter plays a role similar to a SGS model. 

7. naturally adapted to unstructured grid, robust 
enough to cover a wide spectrum of compressible 
flow: from weak linear acoustic waves to strong, 
discontinuous waves (shocks), appropriate for both 
CFD ( computational fluid dynamics) and CAA 
(computational aeroacoustics). 

With an unstructured grid, the CE/SE scheme is easily 
adapted to complicated geometries. More details about 
the unstructured CE/SE method can be found in [13]. 
The weighted a — e CE/SE scheme is used here. 

3 The Impinging Jet Problem 

When a supersonic underexpanded jet impinges nor- 
mally on a flat plate, some shock cell structure is formed, 
and near the flat plate a stand-off or plate shock appears. 
As the jet flow approaches the plate, it is turned into a 
radial wall jet(Fig. 1). Then it is believed that the insta- 
bility waves or vorticities generated in the jet shear layer 
interact with the plate and produce acoustic waves. As 
these acoustic waves propagate upstream to the nozzle 
exit where the shear layer receptivity is the highest, they 
trigger a new cycle of instability waves and thus com- 
plete the feedback loop. 

The impinging jet problems are set up following the 
configuration in Henderson’s experiment [4 - 7]. The 


top boudary 


inflow 


boundary 

outflow 

/ 

boundary: 

nozzle wall 

L \ 




nozzle 
exit 
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sym. axis 


Figure 2: Typical computational domain for impinging 
jet nozzle exit. 

ambient density po, the ambient speed of sound and the 
diameter of the nozzle exit are chosen as the scales for 
density, velocity and length. The computational domain 
includes a near field of the flow. Figure 2 and 3 show 
typical computational domains for the current impinging 
jet problems. As the flow is considered as axisymmetric, 
the two-dimensional computational domains are half of 
the axial section of the corresponding 3-D domains. In 
the case of a small plate with diameter d, equal to D, the 
domain ranges from x = — 3 to x = 7 (sponge/buffer 
zone not included) in the stream direction; while in the 
case of a large plate, the domain ranges from x = — 2 
to x = 3.4 or x = 4.5 with the nozzle exit always lo- 
cated at x = 2. The axial location of the end flat plate 
depends on the ratio of h/D. Here h/D=1.4 and 2.5. Typ- 
ically, there about 35,000 to 40,000 triangulated cells in 
the domain. They are formed by dividing a rectangle cell 
diagonally into four pieces, as shown in Fig. 2 and 3. In 
the area critical to aeroacoustic feedback loop around the 
jet core, the grid sizes are Ax = 0.05 and A r = 0.025 
to ensure enough resolution. 

3.1 Initial Conditions 

Initially, the flow of the entire domain is set at the am- 
bient flow conditions, i. e . , (using nondimensional vari- 
ables) 


NAS A/CR— 2005-2 1 3426 


3 


^Oi""nz 



Figure 3: a: Typical computational domain and grid 
(with buffer/sponge zone at the top) for a large end plate; 
b: an enlargement of the grid around the nozzle exit. 


Pa lj Pa i LI a 0? L) a 0. 

7 

here, the subscript a stands for ’ambient’. 


3.2 Boundary Conditions 

At the inlet boundary, the conservative flow variables 
and their spatial derivatives are specified to be those of 
the ambient flow, except at the nozzle exit, where an 
elevated pressure is imposed, i.e., the jet is underex- 
panded, as in the experiments. As the jet flow at the 
nozzle exit is choked (M e = 1), and the ratio of stag- 
nation (plenum) pressure po to the ambient pressure p a , 
N PR = Po/p a = 4.03, by using the ideal gas isentropic 
relations, it follows that 


Po 

Pe 



1.893, or p e = 2.1289p a 


Other nondimensional flow variables at the nozzle exit, 
with M e = 1 (choked flow with convergent nozzle), are 
given by 

_ 7(7 + 1 )Pe 
Pe 2 T r 

u * = (7^1) > = °> 

where T r is the reservoir (plenum) temperature. We will 
also follow the experimental cold-flow condition where 
the reservoir temperature equals the ambient one, i.e., 
T r = 1. 

At the symmetry axis, i.e., y = 0, a simple reflective 
boundary condition is applied. At the top boundary, the 
Type I CE/SE non-reflecting boundary conditions as de- 
scribed in the next subsection are imposed. The no-slip 
boundary condition is applied on the nozzle walls and the 
end plate. 


3.3 Non-Reflecting Boundary Conditions 

As the spatial derivatives of the conservative flow vari- 
ables are also considered as unknowns, the CE/SE 
scheme supports a simple but robust non-reflecting 
boundary condition (NRBC). Details and proof of the 
new NRBC can be found in [14]. The following is the 
Type I NRBC employed in this paper. 

For a grid node (j, n) lying at the outer radius of the 
domain the non-reflective boundary condition (Type I) 
requires that 


(U X % = (i Uy )] = 0, 

while U™ is kept fixed at the initial steady boundary 
value. At the downstream outflow boundary, the non- 
reflective boundary condition (Type II) requires that 


(U x )? = 0 , 
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Figure 4: Plots 1-8: isobars at different time steps in a 
cycle, showing oscillations of the shock cell, the plate 
shock, the radiating acoustic waves and the vortices in 
the wakes in the small plate case. 

while U "j and ( U y are now defined by simple extrap- 
olation from the nearest interior node j ', i.e., 

u] = u]r 1 (u y )^ = {u y )y\ 

As will be observed later, these NRBCs, when combined 
with the buffer zone, are robust enough to allow a clean 
near field computation without disturbing or distorting 
the flow and acoustic fields. 

4 Numerical Results 

As sketched in Fig. 1, when the jet hits the plate, a 
radial wall jet is formed along with the instability waves. 
It is believed that when the instability waves interact with 
the flat plate and generate acoustic waves that propagate 
upstream via either the jet shear layer or the jet exterior to 
the nozzle lip and complete the acoustic feedback loop. 
In this section, the cases of a small plate and a large plate 
are considered. 

4.1 Jet Impinging on a Small Plate 

Consider a supersonic underexpanded jet from a conver- 
gent nozzle. With a small plate of diameter d placed 
downstream of the jet nozzle exit, the free jet becomes 
an impinging jet. The small plate is indeed a solid cir- 
cular cylinder aligned with the nozzle axis. According 
to the experiment [5], the small plate has a diameter d=D 
and the distance between the jet nozzle exit to the plate is 


Figure 5: Plots 1-8: numerical schlierens showing the 
cyclic movements of the shock cell and plate shock. 

h=1.65D, where D=2.54 cm is the jet diameter at nozzle 
exit. The ratio of stagnation pressure to ambient pressure 
is set at NPR = 4.03. 

4.1 .1 The unsteady oscillating flow Figure 4 
illustrates the isobars at different time steps in about a cy- 
cle, after 1.47 million steps were already run. It is clearly 
displayed that the toroidal instability vortical waves grow 
along the jet shear layer upstream of the plate cylinder, 
and then interact with the bow shock (plate shock) and 
more importantly with the edge of the plate cylinder. The 
vortices are deflected and convect downstream with the 
flow, while the interaction generates acoustic waves that 
propagate in the field. As the waves reach the nozzle exit 
lip, where the receptivity is the highest, another cycle of 
vortices is triggered and the feedback loop is thus com- 
pleted. In our view, the situation is somewhat similar to 
the feedback loop of a high speed cavity noise problem, 
although the situation may be different for a larger plate 
diameter. 

Figure 5 is the numerical schlierens showing the cyclic 
movements of the shock-cell, in particular, the shock cell 
shocks and the plate shock. The strong oscillation of the 
shocks indicates strong, nonlinear waves in the near field. 

4.1.2 Acoustic waves and frequency Fig- 
ure 6 shows the same isobar plots as in Fig. 4 but with 
pressure limited between 0.6 — 0.8 (non-dimensional val- 
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Figure 6: Plots 1-8: isobar contours with pressure limited 
within 0.6-0. 8, showing acoustic waves generation and 
propagation. 

ues). The generation and propagation of acoustic waves 
are clearly observed. The pressure history is recorded for 
a location (—1.5, 2) upstream of the nozzle exit. Fourier 
analysis is then performed to obtain the spectrum. Fig- 
ure 7 is the computed SPL (sound pressure level) plot 
at this location. It is observed that the waves are strong 
(over 150 dB), as expected. The fundamental frequency 
is about 3,300 Hz with ample harmonics. This compares 
well with a similar but not exactly the same case in [4], 
where h=l .65D but NPR=4.40 instead of 4.03, which has 
a fundamental frequency of 3,570 Hz. 

4.2 Jet Impinging on a Large Plate (h/D = 1.4) 

For the case of a jet impinging on a large normal flat 
plate, the plate diameter is d=10D (note that there is 
a buffer zone beyond d=10D in the computational do- 
main), and NPR=4.03. Two spacings between the noz- 
zle exit and the plate, namely, h=1.4D and h=2.5D, are 
considered. Figure 8 is the numerical schlierens at 
different time steps, the shock cell and the plate shock 
(stand-off shock) are clearly displayed but their hydro- 
dynamic movements are much weaker than the previous 
case. Figure 9 illustrates snapshots of a similar experi- 
mental schlierens [6] for qualitative comparison. 

With the pressure values limited in a narrow range of 
0.71 - 0.72, Fig. 10 demonstrates how the acoustic waves 
are generated and propagate in a series of snapshots. Fig- 
ure 1 1 is a snapshot of isobars, giving particular details of 



Figure 8: Plots 1-9: numerical schlierens for large flat 
plate case, the cyclic hydrodynamic movements of the 
shock cell and plate shock are weaker than the previous 
case of small plate. 
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Figure 9: Snapshots of experimental schlierens for a sim- 
ilar case with h/D =1.4, NPR=4.06(from [6 - 7]) 

the vortices along the jet and wall jet shear layers, which 
are believed to be responsible for the noise generation. 
Figure 12 shows the vortices around the jet and wall jet 
in experiment, for qualitative comparison with Fig. 11. 

The computed SPL at a point upstream of the nozzle 
exit is shown in Fig. 13. The spectrum is calculated after 
2.7 million time steps has run. The computed tone SPL is 
slightly greater than 1 10 dB, which means the waves are 
much weaker than in the case with a small plate. Also, 
there are ample harmonics, as described in Henderson’s 
work [6] . 

4.3 Jet Impinging on a Large Plate (h/D = 2.5) 

As the physical behaviors for this large plate case with 
large spacing between nozzle exit and the plate is similar 
to the previous case, only the SPL at a point upstream of 
the nozzle exit is presented in Fig. 14. 

4.4 Comparison to Experimental Acoustic 
Data 

In addition to the above qualitative comparison of the 
structures of shock cell or vortical instability waves 
around the wall jet, comparisons of the computed fre- 
quencies in both small and large plate cases with sim- 
ilar experimental data show that they are similar. The 
comparisons are still somewhat qualitative, since we can 
only find similar experimental data under slightly differ- 
ent operating conditions. The experimental data are from 
Henderson and Powell [4] (figure 1 1 in their paper) and 
Alvi and Iyer [9] (Fig. 22 in their paper). Table 1 shows 
the comparisons. As the data are measured visually from 
the charts, a small error is inevitable. From Table 1, it is 
confirmed that the computed frequencies are in the right 
range. 

Although the sound producing mechanism in the im- 



m c 

C • 0,1 I ♦ 

* o 1 

: 3ji 


# Cl 

IK *' 


— oAa 
l‘% Ol 

* •* 

. ‘<S .. 

1^ Oj 

1 *5$ 





3% <>■ 

. L 

9 t 0< 12 *[2 


(isobars minmax: 0.71-0.72) 


Figure 10: Wave generation in an impinging jet with a 
large plate, h/D=1.4. Isobar values are limited within 
0.71-0.72 in order the weaker acoustic waves can be dis- 
played. 



Figure 11: A snapshot of isobars showing presence of 
various vortical eddies around the jet core and wall jet 
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SPL at (-1 , 2), h=2.5D, NPR=4.03 



Figure 12: Experimental subtracted velocity vectors 
around the jet core and wall jet[6], showing the com- 
plicated vortical eddies there, (note that the flow is in 
opposite direction in the experiment) 


5430 Hz 



Freq. in Hz 

Figure 14: SPL at (-1,2) (upstream of nozzle exit) for a 
large plate, h/D=2.5, NPR=4.03, binwidth = 80 Hz. 


SPL at (-1 ,2), h/D=1 .4, NPR=4.03 



Figure 13: SPL at (-1,2) (upstream of nozzle exit) for a 
large plate, h/D=1.4, NPR=4.03, binwidth = 80 Hz. 


Table 1 : comparison of tone frequencies 



computed 

Henderson [4] 
(exp.) 

Alvi [9] 
(exp.) 

small plate 

3.3 kHz 

3.57 kHz 

— 

d=D 

NPR=4.03 

NPR = 4.40 
S4 mode 

— 

large plate 

4.4 kHz 

5.1 Hz 

4.2 kHz 


h=1.4D 

h=1.4D 

h=1.6D 

(case 1) 

NPR=4.03 

NPR=4.40 

NPR=3.70 
w. lift plate 

large plate 

5.43 kHz 

5.38 kHz 

— 

(case 2) 

NPR=4.03 

h/D=2.5 

NPR=4.40 

h/D=2.5 



pinging jet is not fully understood, the numerical work 
may shed some light on how the acoustic wave is gener- 
ated near a flat plate. Figure 15 shows a typical instan- 
taneous pressure field plus a velocity vector field on top 
of it (Here h=2.5D with a large plate). Near the stagna- 
tion point of the plate, a high pressure bubble with di- 
ameter of about 1-1.5D is formed. Outside the bubble, 
the pressure quickly reduces to the ambient level. In the 
feedback loop, instability waves (vortices) are shed from 
the nozzle lip and grow in strength and size along the jet 
shear layer. When the vortices pass through the tips of 
the shock in the jet core and the plate shock, and enter 
the high pressure stagnation bubble (Fig. 15), the entire 
flow experiences severe changes in terms of all the flow 
variables u, v, p and p. The vortices will be deformed and 
distorted, and an acoustic wave is thus generated. For ex- 
ample, when a vortex turns 90° inside the high pressure 
bubble and exits the bubble along with the radial wall jet 
stream, it interacts with high gradient compression wave 
or shock, and generates an acoustic wave. Outside the 
stagnation bubble, the acoustic wave appears to originate 
from a location near the plate but right outside of the 
stagnation bubble as sketched in Fig. 15. This also ex- 
plains the experimental observation in [6] that the sound 
waves occur at a location of the plate 1.3D from the cen- 
terline (Fig. 12), because this location lies just outside 
the high pressure bubble. 

5 Concluding Remarks 

In this paper, we attempt to numerically simulate the 
important phenomena of a supersonic underexpanded jet 
impinging normally on flat plates, with emphasis on their 
aeroacoustic behaviors. 
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disc shock stagnation 

bubble 

Figure 15: Instantaneous pressure field and velocity vec- 
tor field, h=2.5D. 

The numerical method used is based on the re- 
cently developed CE/SE scheme solving the axisymmet- 
ric Navier-Stokes equations. As the scheme possesses 
low dissipation while being capable of capturing shocks, 
the numerical results compare favorably to both hydro- 
dynamic and acoustic experimental findings [4 - 9] even 
given the two-dimensional axisymmetric approximation. 
However, for jet impinging on a large plate, the sound 
genration mechanism is still not fully understood and 
we have put forward an explanation for the mechanism 
based on the numerical results. More investigations are 
needed in the impinging jet problem to further validate 
the numerical work. 
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